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Abstract 

An universal form of kinetic equation for open systems is considered which naturally unifies classical and quantum cases and 
allows to extend concept of wave function to open quantum systems. Corresponding stochastic Schrodinger equation is derived 
and illustrated by the example of inelastic scattering in quantum conduction channel. 
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stochastic Liouville equation. 



1. Introduction 
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It is seldom when one can say about a quantum system 
Q that it is in a pure state described by certain wave function, 
Q ^{t) . Even closed system generally is in a mixed state and 
should be described by density matrix 1 2 3j, 



> i?(t)= V|*c«W)Pa(*a(t) 
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where "^ait) is a full set of mutually orthogonal solutions 
to system's Schrodinger equation ( (vl>„(t)|v['^(t)) = at 
Q. ^ (3), time- independent Pq, are their statistical weights, 
and whole R{t) is solution to the von Neumann equation. 
All the more if a system S under interest is open, that is 
interacts with some other system W , e.g. "the rest of the 
World" . In such case representation ^ losses meaning and 
density matrix of S becomes "thing in itself" to be found 
|rom von Neumann equation of total closed system S+W . 
' There are many different approaches to this vast problem 
|2I3I4| . However, any of them results in either approximate 
kinetic and stochastic equations for S or formally exact but 
too complicated ones. In present paper, in part following 
|5|6|7|8|9^0] . one more approach is described which leads 
to such kinetic equation which is at once simple and exact. 
Its other advantages are intuitive obviousness, explicit sep- 
aration of contributions from S and W and possibility to 
reformulate it in terms of wave function of S itself, thus re- 
habilitating representation ([T]) as applied to open systems. 



Email address: kuzovlevSkinetic.ac.donetsk.ua 

(Yuriy E. Kuzovlev). 



2. Stochastic treatment of dynamic interactions 

As usually, let us divide Hamiltonian H of the whole 
system S+W into three parts and assume that the inter- 
action part. Hint , has bilinear form: 



H — Hs + H\Y + Hin 



(2) 



where operators Sj and Wj act at different Hilbert spaces 
of S and W , respectively. Such a decomposition of the in- 
teraction Hamiltonian is always possible (and, as a rule, it 
directly follows from physics of interaction) . For any oper- 
ator O define Liouville and Jordan super-operators: 

£{0)A = ^{OA- AO) , I1{0)A = ^{OA + AO) , 

where A is arbitrary operator. Thus C{H) is the Liouville 
evolution operator of S+W . One can see that correspond- 
ingly to ^ it decomposes to 

C{H) = C{Hs + Hw + E SjWj ) = (3) 

= CiHs) + C{Hw) + E ['^iS,) n{W,) + Tl{S,) C{W,) ] 

From the viewpoint of S it is formally equivalent to 

Cs{t) ^ CiHs) + J2[x,it)C{S,) + y,{t)n{S,)] (4) 

where Xj (t) and i/j (t) play as arbitrary scalar time func- 
tions although in respect to W they are super-operators 



X-j (t) 



_ -C{Hw)t 



C(Hw)t 
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In detail, if p{t) denotes joint density matrix of 5+W , 
which undergoes the von Neumann equation p — C{H)p, 
and ps{t) = Trvi/p(i) density matrix of 5, then 



Ps{t) = Tr^/ 6xp 



I 

j Cs(t')dt' 



P(0) 



(5) 



where 6xp means chronologically ordered exponential. 

These facts prompt us to consider Xj{t) and yj{t) as 
random processes while operation Tr w in ([5]) as statistical 
averaging over them. In order to realize such "stochastic 
representation of dynamic interaction" between S and W 
in a simple unambiguous way, let us assume that at some 
initial time moment , e.g. to = Q , before the interaction 
S and W were statistically independent one on another: 
p{io) — Ps^Vm^^ ■ Then we can write 

ps{t) = {R{t)) , (6) 
where R{t) satisfies "stochastic Lioville equation" 



dR{t) 
dt 



Cs{t)R{t) 



(7) 



with initial condition R{to) = Pg"^ , and angle brackets 
(...) designate the averaging. At that, evidently, for any 
statistical moment (...) —Ttw ■■■ Pyy\ where 
symbolizes chronological ordering, and complete statistics 
of random processes Xj(t) and yj(t) is summarized by 
their characteristic functional as follows. 



exp y ^ [ u, (t) X, (t) + f, (t) y, it) ] dt ) 



Cwit) dt 



(in) 

Pw 



(8) 



Cwit) ^ CiHw) + Y.[uj{t)Il{Wj) + fjit)CiW,)] 



Here Ujit) and fjit) are arbitrary probe functions, the 
integrations begin at t — tQ , and, if necessary, tQ can be 
moved away to — cx) . 

Physical meaning of Xj it) in (g]) and ©-([T]) is quite clear 
from (|3]): Xjit) merely reproduce Hamiltonian perturba- 
tion of 5 by W as prescribed by ([2]) . In order to understand 
meaning of yjit) there, notice that according to ([8]) the 
super-operators n(Wj) represent observation of variables 
Wj in subsystem W by subsystem S . Hence, analogously 
and symmetrically, super-operators n(S'j) in Csit) cor- 
respond to observation of variables Sj in iS by W while 
yjit) , according to ([8]), represent conjugated Hamiltonian 
perturbation of W by 5 . Therefore yj it) introduce also 
response of W to this perturbation and eventually various 
feedback effects, in particular, dissipation and renormaliza- 
tion of dynamic properties of S . 

By these reasons, generally Xjit) in ([7]) behave like 
conventional random processes, but yjit) are extremely 
unusual ones. Indeed, if Ujit) = in ([8]) then Cwit) turns 



into purely Liouville super-operator, therefore the char- 
acteristic functional turns into unit regardless of fj it) . 
Thus any statistical moment of yjit) themselves equals 
to zero, {yiti)...yitn)) = . Moreover, any moment where 
some yjit) has most late time argument also is zero: 
(a;(Ti)...a;(T„)y(ii)...y(i„)) = if maxj(tj) > maxj(Tj-). 
However, in general (a;(r)y(i < r) ) ^ 0. Physically, 
these features of yjit) trivially mean that perturbation 
of W by iS (and conjugated observation of iS by W ) al- 
ways foregoes response of W to it, in agreement with the 
causality principle. At that, (a;(T)?/(i < r) ) and other 
cross-correlations just represent dissipation in S . 

In the classical limit super-operators Ci-..) and H(...) 
become the Poisson bracket and mere multiplication by a 
phase function, respectively, but physical sense of and 
yit) and their main statistical properties do not change. 
More details about them and various generalizations of the 
"stochastic representation" can be found in |5I7I8I10| . 

3. Stochastic Schrodinger equation 

If initially (before the interaction) S was in a pure state, 
pf''> = |$(™))($(™)| ^ then equation ^ has a solution 
which all the time keeps such form: i?(i) = |^'(i))(5'(t)| , 
where ^'(i) satisfies stochastic Schrodinger equation 



with initial condition 5* (to) = ^ , and 



(9) 



Wjit) = Xjit) + ihyjit)/2 

are mentioned as complex- valued random processes. Hence, 
solution to ([7]) under arbitrary initial condition i?(to) , let 
^{m) _ I^q"^) Pa (*q"^| in the diagonal form, can 
be written in the form ([1]) where each of ^ait) evolves 
according to ^ but all Pa stay constant. 

Because of complexity of Wj it) , the evolution governed 
by ([9]) is not unitary: 



dt 



(*a(t) 1*^(0) - (*c.(t)|^,|*M^)) ^ (10) 



Nevertheless it is unitary on average: 



d 
dt 



{^^it)\^f,it)) ) = {J2 %(t)(*a(t)|5,|vl/^(t)) ) = , 
( (*a(t)|*0(t)) ) = i^Ml^pito)) « Sa(3 (11) 



This statement directly follows from the above mentioned 
specific statistical properties of yj it) . 

Let the interaction by its nature is time-localized, that 
is has a character of scattering. In such situation for any 
particular "stochastic pure state" one can write 

\^it)) = 1*0(0) + l^^it)) 
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with [^"(t)) describing free evolution of S and |^'*(t)) at 
t ^ oo representing result of the scattering. Then due to 
(dlD ( ) = (*"|*°) and consequently 



(12) 



This is nothing but "the optical theorem on average" which 
equally holds for both elastic and inelastic scattering. It 
demonstrates also that average wave function can 
bring useful information about an open quantum system. 

4. Interaction with thermal bath 

Of special interest are cases when W is very large system 
in thermodynamical equilibrium, p'^^ cx exp(— iJ^/T) , 
so that all xj (t) , i/j (t) and Wj (t) are stationary random 
processes. If one assumes (without loss of generality) that 
Ttw Wj py^^ = , then all the processes have zero mean 
values and all their pair correlations reduce to function 

Kjm{T) = TrvF Wj e:xjp{-iTHw/fi)Wm p'^^ 
One can routinely verify that 

{w,{T)wmm = {w*{r)w;Mr = 

= K,ra{r)0{T)+K*„^{T)6{-T) , (13) 

ifJ,^(T) ^ {x,{T)xrM) - Reif,„,(r) , 

where 6{t) is the Heaviside step function. The last equality 
expresses equilibrium of W and in its turn implies definite 
relation between K^^{t) and Kj^{t) : 



Kj,niT) = J cos{ujt) ajm{uj)duj , 





(14) 



K'^y (t) 



2e{T) 



sin(wT) tanh 



2T 



where ajmiw) is non- negatively defined spectrum matrix. 
In essence, this is usual fluctuation-dissipation theorem. 

Of course, the W 's equilibrium results also in definite re- 
lations between higher-order correlators of x{t) and y{t) . 
They were discussed in [11]. In the classical limit all they 
can be summarized by symbolic equality 

, , , , 1 dxj(t) 

e,y,{-t) X y,{t) + if , 

where x means statistical equivalence under time inver- 
sion, t ^ —t , and Ej ± 1 indicates time parity of variable 



Xj{t) from viewpoint of Hamiltonian mechanics, so that 
Xj{—t) X £jXj{t) . Correspondingly, the second-order rela- 
tion expressed by (I14p in the classical limit simplifies to 



K'-'y (r) 



T 



il A. ft-^^ 

dT ^'^ 



(r) 



5. Inelastic scattering center in conducting channel 

Let iS be one-dimensional conduction channel formed by 
discrete sites (n=... — 1,0,1,...), W be thermostat, and 
their interaction realizes through randomly varying poten- 
tial localized at only one site, for certainty, n = . Then 
{Hs)n'n = e [2(S„'„ - Sn'n-1 ~ Sn',i+i ]/2 , where e is chan- 
nel's bandwidth, and '^Wj{t)Sj — w{t) S , where 5'„'„ = 
'^no'^ii'o ■ We want to consider evolution of wave function 
5'(<, n) which initially, ai t ~ to < , represents a wide 
wave packet placed on the left from n = and moving to 
the right with momentum hko > . Thus let us divide it 
into free and scattered parts: ^(i, n) — '^^{t, n) -f^''*(t, n) , 
where vi/° = - iHs ^"/h , so that jH) transforms to 



d'^''it,n) 
Jt 



w{t)dno'^it,0)] , (15) 



with initial condition ^'^(toi n) — . 

Strictly speaking, such division presumes that particle 
scattered by fluctuating potential can not be captured by 
it. If it is really so, that is ^''(t, rt) ^ at i — > oo for any 
given flnite n (first of all, n = 0), then at large enough 
t > after scattering we can determine probabilities of 
reflection, 5ft , and transmission, T , as 



r 



n< 



1*1 



n> 



s|2 



n< 



2 ^ 



(16) 



Here we took into account mirror symmetry of scattered 
wave, ^''*(t, — n) = ^''^(i, n) , which obviously follows from 
. After averaging over the thermostat and combin- 
ing with general "optical theorem" (fT2|) we have 

(5R) = 1-(T) = -ReE*°*(^,")(*'(^,")) (17) 

n 

Alternatively, we can describe the scattering in terms of 
outgoing wide wave packets with discrete wave numbers fc 
(separated by a width of ingoing packet in reciprocal space) 
and amplitudes = ^\ + *| . Then = , and 
optical theorem (fT7|) takes form 



(5R) 



-Re ) 



(18) 



The second equality here is again consequence of the mirror 
symmetry of scattering which implies *1 y. = 4'^ . 

Next, let us divide 5R and T into elastic and inelastic 
parts marked by "el" and "in" , respectively, so that 



(19) 
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and use obvious general inequalities 



6. Conclusion 



By applying them to (|19[) and combining with (jl8[) it is 
easy to obtain 



(T*"} = (T) - (T'^') < (T) - (r)2 



(20) 



In fact these are identical inequalities since = T'" 
already due to the mirror symmetry. Their sum gives us 
restriction on total probability of inelastic scattering, V™ : 



(K™} + (r™} < 2 (SR) (1 - (SR)) < 1/2 (21) 



This seems interesting and important consequence from the 
"unitarity on average" . 

More careful consideration bases on stochastic integral 
equation for ip{t) = ^(t,0) which follows from ([TS]): 



(22) 



i^it) = ^"(t,0)-- / Go{t^t',0)w{t')^{t')dt' , 



to 



Go(t, n) — / exp[mfc— iE{k)T/h] 



dk 
2^ ' 



with E{k) = e [ 1 — cos k]/2 being dispersion of the chan- 
nel. Let us confine ourselves by weak scattering, in the 
sense that (w(t)) = and (7{ll!) <^ he , where spectrum 
cr(w) of the fluctuating potential is mentioned like in (fT4| 
as spectrum of its real part x(t) = Re u'(t): ( a;(T), a;(0)) = 
cos(tJT) <7{uj) duj . Thus, according to P^ - (f]7 



{w{t)w{0)) = K{\t\) , {w*{t)w{0)) = K{t) , 



K{t) 



e'"^ +exp(Sw/T)e-'''^^ 
l + exp{huj/T) 



a{uj) doj 



(23) 



Then simple analysis of equation l|22|) in the correspond- 
ing standard "one-loop" approximation (or, in other words, 
"the Bourret approximation" [12]) yields 



C 



l + C 

-Eo}/h 



C 



1 



<y{\i^\) dio 



(24) 



-Eo/h 



[1 + exp (hoj/T)] ^ [Eq + nw){e - Eq - hw) 



where Eq = E^ko) . In fact, due to our assumption, C <^ 1 
and at that, naturally, almost all the scattering is inelastic, 
except may be case of very small particle's energy Eq (for 
more details see [9]). Contributions to the integral from 
Lo > and w < correspond to scattering with energy 
radiation or absorption by the thermostat, respectively. 



We just had a chance to make sure that there exists uni- 
versal form of stochastic Liouville and Schrodinger equa- 
tions which is exactly valid for arbitrary open quantum 
system S representable as a part of some closed system 
S+yV . Simultaneously there exist universal exact recipes 
for determination of statistical characteristics of random 
sources, or noises, w{t) , what enter these equations. 

Peculiarity of thus stated "stochastic representation of 
dynamic interactions" (between S and W ) is that stochas- 
tic images, w{t) , of Hermitian operator variables W{t) of 
W behave like a c-number (purely commutative) complex- 
valued random processes. Their imaginary parts, y{t) = 
Im w(t) , introduce dissipation into 5's evolution by means 
of virtual violation of its unitarity. Nevertheless, due to 
specific statistical properties of y{t) , statistical averaging 
restores the unitarity although keeps the dissipation. 

It is useful to underline that our approach to open sys- 
tems makes it possible to obtain final results under inter- 
est without conventionally thought "kinetic equations" , i.e. 
omit stages of approximate derivation of a kinetic equa- 
tion and then its approximate solving. Illustrations were 
done in [516171819] and just above. From the other hand, our 
approach gives all tools for correct construction of kinetic 
equations as well as "Langevin equations" for 5's variables 

m- 

Importantly, Hamiltonian character of joint dynamics 
of S+W implies general statistical connections between 
real and imaginary components of w{t) (especially when 
W is a thermostat [H]), which help to convert simple ab- 
stract models of stochastic processes into our approach. At 
that, one can exploit all the theory of linear dynamic sys- 
tems with randomly varying parameters (including theory 
of waves in random media |12j ) . 
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